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ABSTRACT 

A Poisson equation on a rectangular domain is solved by coupling two 
methods: the domain is divided in two squares; a finite element approximation 

is used on the first square, and a spectral discretization is used on the 
second one. TVo kinds of matching conditions on the interface are presented 
and compared; in both cases, error estimates are proved. 
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1. INTRODUCTION 


To approximate the solutions of partial differential equations, a number 
of methods can be successfully applied: among them, spectral type methods, in 
which the discrete solution is a polynomial of high degree, are known to be 
very accurate when the solution approximated is very smooth (see [GO][CHQZ] 
for a general description of these methods). Their main drawback lies in the 
difficulty to take into account the singularities of the function approximat- 
ed, as well as the difficulty in handling domains with a complicated bound- 
ary. This last problem is usually solved by decomposition into subdomains 
and/or transformation of coordinates. On the other hand, the finite element 
method, where the discrete solution, restricted to very small domains called 
"elements,” is a polynomial of low degree, is well suited to problems with 
complex geometries, but its accuracy is limited by the degree of the polynomi- 
als (general properties of finite elements are analyzed in [C]). Several at- 
tempts have been made to set the two methods in a unified framework and obtain 
the advantages of each one. The spectral element method [P], which consists 
of using a spectral algorithm on a fixed number of subdomains, is presently 
developed for a growing number of problems (see, for instance, [F], [KP], and 
[MP ] ) ; on the opposite side, the so-called p-version of the finite element 
method, where the discrete functions are polynomials of fixed high degree on 
each element, is studied by several authors ([BSK], [SV] , [V] for instance). 

The idea of this paper is very different: as previously presented by 
K. Z. Korszak and A. T. Patera [KP] , it consists of dividing the domain where 
the problem must be solved in two parts; then, the problem will be approxi- 
mated by a finite element method on the first part and by a spectral method on 
the second. Consequently, the discrete space will consist of functions which 
are piecewise polynomial on one part and a restriction of a high degree poly- 
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noroial on the other with a matching condition on the interface. Here, we pre- 
sent and compare two kinds of matching conditions: the first kind is a ponc- 

tual one, i*e., we bias the functions to be continuous at the nodes of the 
finite elements on the interface; the second kind is an integral one, since we 
require the trace of the finite element function on the interface to be the 

o 

L^-proj ection of the trace of the polynomial onto the finite element space. 
Of course, both algorithms will be nonconforming in the general case, since it 
is impossible to match a high-degree polynomial and a piecewise polynomial 
function on the interface in a continuous way. However, in a finite element 
context, nonconforming methods have proved themselves to be as efficient as 
the conforming ones (see for instance [CR] or [RT]). Moreover, numerical ex- 
periments* [KP] already show the interest of the coupling technique, which 
turns out to be easy to implement and very flexible to fit both the problem 
and the domain. 

In this paper, we analyze the coupling method on a test problem and in a 
model domain. The domain is simply the rectangle ft - (-1 , l)x(0, 1) , which 
we divide in two parts, ft” = (-l,0)x(0,l) and ft + = (0,l)x(0,l); we de- 
note by Y the interface {0}x(0,l), and by n the unit vector which is 
orthogonal to Y and directed from ft” to ft + . For a given function 
f on ft, the Poisson problem we want to approximate is the following one: 
Find a function u on ft such that 


( 1 . 1 ) 


-Au = f in ft, 
u 55 0 on 3ft. 


1 N. Debit, 


Thesis in preparation. 


n 


Figure 1.1. Decomposition of the domain. 

An outline of the paper is as follows. In Section II, we introduce the 
discrete spaces and state the discrete problems. Sections III and IV are de- 
voted respectively to the analysis of the consistency error and of the approx- 
imation error. The final error estimates, together with concluding remarks, 
are given in Section V. 

The main results of this paper were first presented in [BDM]. 

Notation: Let A denote any open interval of the real line or any domain in 
with a polygonal boundary. For any real number s, we consider the 
classical Hilbert Sobolev spaces H S (A), provided with the usual norm 
H.B g A , and also, when s is an integer, with the semi-norm |.| g For 

any real number s ^ 0 and any p, 1 _< p +°°, we also use the Sobolev 

spaces W S,P (A), provided with the norm J.l g Finally, for any real 
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number s ^ 0, H®(A) stands for the closure in H S (A) of the space of 
indefinitely differentiable functions with a compact support in A. 

Throughout this paper, with any function v defined on ft, we associ- 
ate the pair v = (v ,v ), where v (resp, w) denotes the restriction of 
v to ft (resp. ft ). The following scalar product on L (ft )xL (ft ) 

(1.2) (u ,v ) = / u (x)v (x)dx + / u + (x)v + (x)dx 

ft - n 

2 

coincides with the usual one on L (ft) . We also provide the product 
H* (ft - )xH*(& + ) with the norm 

(1.3) II v II = [(v ,v ) + (Vv ,Vv )] 1/z ; 

the space of pairs v in H (ft )xH (ft ) with v continuous through y> 
is isomorphic to H^(ft). Finally, we define on H*(ft~)xH*(ft + ) the bi~ 

linear form 

(1.4) V (u*,v*) e [H 1 (ft _ )xH 1 (ft + )] 2 , a(u*,v*) = (Vu*,Vv*). 

2 

Clearly, for any function f in L (ft) , problem (1.1) is equivalent to the 
following one: Find u in H^(ft) such that 

(1.5) V v e H^(ft), a(u*,v*) = (f*,v*). 

This variational form is precisely the one which will be used in order to 
define the discrete problems. 
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In all that follows, c, c ' , c" are generic positive constants, 

independent of the discretization parameters. 


2. THE DISCRETE SPACES AMD PROBLEMS 


2.1. Definition of the discrete spaces . 

We have to define a discrete space on each subdomain ft and ft + , 
and then we must match conditions on the interface. 

Let h be a real parameter, 0 < h £ 1, which will tend to 0. With 

each value of h, we associate a triangulation of the domain ft , i.e., 

a finite set of triangles such that the intersection of two triangles is 
either empty or a vertex or an edge and such that 

(2.1) TT = K; 

h 

h is the upper bound of the diameters of the triangles of We denote 

by h^ the diameter of any triangle K in T^, and by the diameter of 

the inscribed circle in K. Next, we assume that the family ( is 

regular in the following sense (cf. [C, Section 3.1] or [B, Def. 3.1]: there 
exists a constant t > 0 such that, for any h, and for any K in 7^, the 
following inequality holds 


( 2 . 2 ) 


p k 1 T V 


- 6 - 


Let k be a fixed integer _> 1. For any closed subset A of It 
( resp . B?), we denote by P^A) the set t ^ e restriction to A of 

polynomials of one variable (resp. two variables) with total degree k. 
With any triangulation 7^, we have the associated finite dimensional ^ 
defined by 


(2.3) ^ = (v h e C°(fi ); VK e T h , v h | R e * k (K) and v h = 0 on \y}. 
We also need the finite dimensional trace space 


(2.4) x h = {v h j Y , v h e x h >. 

In order to build an appropriate basis of and we consider each 

triangle K as the support of a Lagrange finite element (K, P^OO,:^), 

where E is the set of all points in K with barycentric coordinates 

K. 

i/k, j/k, and (k-i-j)/k, 0 £ i , j k, i + j k; it is well-known [C, Thm. 

2.2.1] that this set of points is P k (K)-unisolvent . Next, we set 

< 2,5) \ = MceT. a K’ 

n 


and al so 

(2.6) 5 h = (a e E h Hy}. 

To each point a in O (£TUy) , we associate its Lagrange interpolant, 

i.e., the unique q a of which is equal to 1 in a and vanishes at any 


other point of 5, . Then, the set {q ,a e (_Jy)} is a basis of 

n an 

the space and the set (q a |^.,a e is a basis of the space x^. 

Next, let N be an integer >1 1* which will tend to +". For any 

integer n ^ 0, we denote by Q n Cft^ ) the set of the restrictions to 7^ 
of polynomials of two variables with degree n with respect to each vari- 
able. For each integer N, we consider the finite dimensional space 
defined by 

(2.7) ^ = {v N e (^(G*); v N = 0 on 30 + \y}. 

Let (L ) __ be the family of Legendre polynomials on [0,1], i.e., of or- 
thogonal polynomials on [0,1], such that 1^, n e 1) Is of degree n and 
satisfies L n (0) = 1. We recall that the set { L m © L n > 0 < m,n < N} is a 
basis of Q n (H + ). However, we shall also characterize the polynomials of 

S N by ponctual values: let 0 _< j _< N, be the roots of the polynomi- 
al C ( 1 — C)L n (C), with 0 = < ... < ? N = 1; we set 

(2.8) S N = {(? i ,C j ), 0 < i,j < N), 
and 

5 n = S n Hy - {(O.Cj), 1 < j < N-l } . 


(2.9) 


Figure 2.2. The triangulation and the set 


S N* 

Finally, with each value of h and N, we associate the discretization 

"1 

parameter 6 = (h,N ). The pair u , where u is the solution of problem 

(1.1), will be approximated in a subspace of made of a11 pairs which 

satisfy a matching condition on the interface y. More precisely, we are 
going to consider two kinds of matching conditions, with which we associate 
two kinds of discrete spaces, both denoted by V^: 

1) ponctual matching condition: the space is defined by 


( 2 . 10 ) 


V 6 = {v 6 " ( VV £ *h X V Va e «h» v h (a) = V a)}; 
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2) integral matching condition: the space is defined by 

(2.11) V* = (v* = (v h ,v N ) e ^ x Vq h e x h , / (v h - v N )(0,y)q h (y)dy = 0}. 

p 

The space has already been used in [KP], Here, we compare the two 

kinds of spaces, and the integral matching condition will turn out to be 
better. 

Remark 2. 1 : We immediately notice that both methods are nonconforming since 

* 

the functions v^ , associated with a pair of are generally 

discontinuous through y and consequently do not belong to Hq(J 2). In- 
deed, for N _> k, a function v^, with the pair v^ in X^xX^, 
belongs to Hq(£ 1) if and on ^y if it is a polynomial of P^(Jl)OHQ(ft) . 

Remark 2.2: From a numerical point of view, to enforce the ponctual matching 

condition, one has to interpolate polynomials of at every point of 5^ 

and hence must store the values L (a), 0 < n < N, a e On the other 

n — — n 

hand, to enforce the integral matching condition, one needs to store the 

integrals / L (y)q (0,y)dy, 0 < n < N, a e £, . Consequently, the cost of 
y n a — — h 

the tvro methods is of the same order. 

However, when k is equal to 1, for a given value of N, it is possible 
to choose the triangulation such that the sets and £ N 

coincide. Then, since the polynomials of Xjj are characterized by their 
values at the points of S^, enforcing the ponctual matching condition would 
be less expensive. But this would require very strong restrictions on the 
triangulation T^; in particular, the parameters h and N would be linked 
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by a relation of the type h >_ cN -1 . Moreover, the triangulation could not 
be uniformly regular since it is well-known (see [S, Thm. 6. 21. 3]) that the 
c j» 1 < i < N-l, satisfy: ^ = sin 2 0^., with (2j - l)ir/4N < 9^ < (j + 

1)it/ 2(N + 1); hence the points of are not all equally distributed 

(they cluster to ±1). That is why we would not recommend such a choice. 


2.2. Definition of the discrete problems . 

We are now in a position to define the discrete problems. We recall (see 
[DR, Section 2.7] or [H, Chapter 25]) that there exist positive weights 
p , 0 < j < N, such that the Gauss-Lobatto quadrature formula 


( 2 . 12 ) 


1 N 

/ «(5)dC « l $(C.)P, 

0 j=0 J J 


is exact on all polynomials with degree _< 2N - 1. 

With each point a = (? ^ ,C j ) in S N , we associate the weight 
P = P.P,« We now introduce the following discrete bilinear form on 

a i 3 
L 2 (fi~)xC°(£f ) 

(2.13) (u*,v*) 6 = I KeI J K u"(x)v"(x)dx + £ ae5 u + (a)v + (a)p a , 

h “N 


which coincides with the usual scalar product on l? (ft ) x <^j_jCfl f )» Finally, 
we define on H*(ft )xC^(?2 + ) the bilinear form 


V (u*,v*) e [H 1 (fi”)xC 1 (fi + )] 2 , a.(u*,v*) = ( Vu * , Vv*) . . 


(2.14) 
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Then, for any pair f* given in for each kind of match- 

* 

ing condition, the discrete problem is the following: Find , with 

in Vg, such that 

( 2 * 1 5 ) V v* e ¥ a , = (f*,v*) 6 - 

Remark 2.3: Of course, in definition (2.13), one could, by using a quadrature 

formula, replace each integral J^u - ( x) v - ( x) dx by its approximation. The 

o 

resulting algorithm will be thoroughly analyzed. 

We recall [CQ, Lemma 3.2] the property 

N ? _l 1 9 

(2.16) l L (e > p = (2 +N l ) f L (C )dC. 

j=0 3 J 0 

Since the quadrature formula (2.12) is exact on all polynomials of degree 
_< 2N-1, the discrete scalar product (.,.)^ is uniformly equivalent to 
(.,.) on l (£1 - )xQ^(TF ). Consequently, the form satisfies the 

following properties of continuity 

(2.17) V (u*,v*> e [H 1 (fi _ )x<^(lf )] 2 , |a 6 (u*,v*)| < cHu*H Hv*H , 

and of ellipticity 


^N. Debit, Thesis in preparation. 
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(2.18) Vue H 1 (a")x<^(n ), a 6 (u ,u ) > (Vu ,Vu ). 

Now, since both 3ft - n9ft and 3ft + n9fl have a positive measure, it 
follows from the Poincare-Friedrichs inequality that the semi-norm: 

* k k 1 /? 

v -*■ (Vv ,Vv ) is a norm equivalent to II . II on the space 

{v*eH 1 (Jl + ); v = 0 on 3ft} , which yields in our particular case 

(2.19) V u* e (ft~)xQ^ ) , a fi (u ,u ) _> clu » 2 . 

Thus we have proved the following result. 

Proposition 2.4 : In both cases of ponctual and integral matching conditions, 

k 

problem (2.15) has a unique solution u^ with u^ in V^. 

The purpose of what follows is to give an error estimate between the 
solutions u and u fi of problems (1.1) and (2.15) respectively. We begin 
with a classical bound (see [C, Thm. 4.2.2]). 


Proposition 2.5: In both cases of ponctual or integral matching conditions, 
the solutions u and u^ of problems (1.1) and (2.15) satisfy 


k 

llu 


- Ug II £ c inf ^ {llu - Vg If 
v 6 eV 6 


k k k k k 

+ sup . [a(v.,w.) - a.(v *,w.)]/ Iw.l} 


w 6 eV <5 


* * 


k k 


+ sup * [(f ,w 5 ) - (f ,w (S ) (S ]/lw (S ll 


w 6 eV 6 


( 2 . 20 ) 



- 13 - 


+ sup A 

w 6 = (w h> W N )eV 6 y 


/ ( 3u/9n) (0,y)(w N - w h ) (0, y)dy/ llw^ I 


7C 

Proof: Let v g be any element in V fi . Using (2.19), we have 


* * 2 * * * * 

C »V V < a 6 (v V U 6' V 


* * * 


* * * 


* * * 


* * 


= -a(v 6 ,u 6 - v 6 ) + a(v fi ,u 6 - v fi ) - a^v^- v fi ) + (f ,u fi - v fi ) fi . 


Next, it follows from (2.1) that, for any w* in H 1 (£T)xH 1 (£2 + ) such 
that w vanishes on 3ft, 


(f*,w*) = / f(x)w(x)dx = -/_ (Au)(x)w (x)dx - / + (Au) (x)w + (x)dx 
£1 ft ” £1 

= a(u*,w*) + / (3u/3n)(0,y)(w + - w )(0,y)dy. 
Y 


^ 

Setting w = u^ - v^ = ( w h » w N ) and combining the result with the previous 
inequality, we obtain 


cllu 6 


*2 , * * * * , * * * , * * 
V < a(u - v 6 ,u 6 - v 6 ) + a(v 6 ,u 6 - v fi ) - a 6 (v 6 ,u 6 




* 

* u 6 


- v 6 ) + 


(f ,u 



+ / 
Y 


(3u/3n)(0,y)(w N - w h )(0,y)dy, 


and (2.20) follows. 


We are now interested in deriving a bound for: 


-14- 


1) the consistency error term 

sup * / (8u/9n)(0,y)(w N - w h > (0,y)dy/ llw* II , and 

w 6 = (w h’ w N )eV 6 Y 

2) the approximation error term 

ie * -k 

W u -v 6 11 

since estimating the two other terms of (2.20) is a standard result in 
spectral methods. 


3. ANALYSIS OF THE CONSISTENCY ERROR 

The aim of this section is to study the term 

/ (3u/9n)(0,y)(w N - w h )(0,y)dy, for any pair = (w h ,w N ) in V fi , 

Y 

where u is a given function on 0 which we shall assume to be 
sufficiently smooth. This analysis involves only one-dimensional approxima- 
tion operators. 

3. 1 The case of ponctual matching condition . 

We recall that there exists an interpolation operator from 

{v e C° (y) ; v(0) = v(l) = 0} into such that, for any function v 

continuous on T and vanishing at 0 and 1, is the only element of 

Xft which satisfies 


(3.1) 


V a e (c h v)(a) = v(a). 
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Mo re over , for any real number £, l_<£_<k+l, there exists a constant c 

such that, if the function v belongs to W £,p (y) for a number 

p, 1 < p jC +°°, the following interpolation error holds [C, Thm. 3.1.5] 


(3.2) 


Sv - -t. vfl 

h m,p,Y 


< ch £-m lvl . , m = 0 or 1, 

- ^,p,Y’ 


First, we prove suitable inverse inequalities in the space <^(0^), 

which complete those given in [Q, Section 3.1] and could be useful in other 
applications . 


Lemma 3.1: Any polynomial on [0,1] with degree N satisfies 


(3.3) 


"%*!,!, (0,1) - 2 N "Vo,«,(0,l)* 


Proof: First, we have 

"Vo, i, (0,1) - "% "o,°°, (0, 1 ) ’ 

Next, let a x be the zeros of q' (i.e., the extrema of q N ) 

belonging to (0,1), in increasing order. Of course, K is < N. Then, it 
suffices to compute 

/ iy°! d? 

0 

K-2 

- w - y°)i + i i q N (< w ■ v\ ) i + i q N (l) ■ Vvi’i 

k— 1 


< 2 KI| q N n 0> „ >(0}1) * 
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Remark 3.2 : The inequality (3.3) is optimal. Indeed, choosing q^ equal to 

the Chebyshev polynomial with degree N, i.e., T N (c) = C0S ( N Arcos(l - 2^)), 
we see that 

1 

/ tooldc = 2N . 

0 

Corollary 3.3 : There exists a constant c such that, for any p and p , 

1 <_ p 2 and 1/p + 1/p* * 1, any polynomial q N on [0,1] with degree 
< N satisfies 


(3.4) 


lq.,1, < cN 3 2/p lq„l , 

» 1,P,(0,1) - N 0,p ,(0,1) 


Proof: We recall the well-known inverse inequality [CQ1, Lemma 2.1], 


" q N H 1, (0 , 1 ) — CN “Vo, (0,1)* 

Interpolating between this inequality and (3.3) gives (3.4), with a constant 
independent of p (see [BL, Thm. 1.1.1] or [LM, Chap. 1, Th . 5.1]). 


Corollary 3.4 ; There exists a constant c such that, for any real number p, 
1 < p < 2, any polynomial qjj on [0,1] with degree < N satisfies 


(3.5) 


" q N"l,p,(0,l) i cN / /P|q N l 1/2 , (0 , 1 ) 


Proof: For p > 1, setting 1/p* = 1 “ 1 /p , we know from [BL, Thm. 6.4.5] 

that H^^(y) is the interpolation space with index 1/2 between W^ , 1 > (y) 
and L p *(y); hence we derive (3.5) for p > 1 from (3.4). That implies 
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%\i,(0,i)< c ^ /2 ' 1/P, Vi/2,(0,i> 

for p tending to 1, which proves (3.5) with p - 1. 

We need a precise version of the Sobolev imbedding (further details can 
be found in [T] for instance) . 

Lemma 3.5: There exists a constant c such that, for any real number p, 

2 < p < +°», any function v in H^^(0,1) satisfies 

(3 ’ 6) ,v, 0,p,(0,l)i c/?,v, l/2 > (0,l)- 

Proof: Since H 1 / 2 (0,1) coincides with the space of the restrictions of the 

functions of H^ 2 (3R), it suffices to prove (3.6) with (0,1) replaced by R 
First, let e be a real number > 0. For any function v in H*'^ +e (j*), 

A 

denoting by v the Fourier transform of v, we have 
v(t) = (2ir) -1 ^ / v(Oe -i ^d5 , so that 

,v, o,- > i t < c 4 i ;(c) ' dE 

< c(/ |v(C)| 2 (l + |5| 2 ) I/2+5 dC) 1/2 (J (1 + |C| 2 r 1/2 ‘ E d5) l/2 . 

m m 

Since the map: v (/ |v(C)|^(l + 1 5 1 ^) ^^ +e d5 )^ ^ is a norm on (H) 

equivalent to the one which is obtained by interpolation, we must estimate the 
constant / (1 + |C | 2 ) -1 ^ 2-e d£. The change of variables 5 = cotanm 


gives 
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f ,1 .. | _ | 2. — 1. /2— e . -l+2e . / /2 -l+2e , 

J (1 + |5| ) d£ = 2 J (sina)) dw < c J a) dw=c/e. 

m oo 


We deduce 


|V, 0,-,]R < c/ ^ ),v, l/2 +£llt 


Finally, for a real number p > 2, we choose 5 = l/(p - 2) and we interpo- 
late this last inequality with the identity llvB A _ = llvll A . We obtain 

U , JK U , K 


*0,p,H < <c//F)1/<I+2e>,v, l/2,B < ,v, l /2, ]R ' 


We are now in a position to prove the following result. 


Proposition 3.6 : For any function u in Hq(A) O H 2 (ft) , the following 

estimate holds for any w£ = in 


(3.7) / (3u/3n) (0,y)(w M - w. )(0,y)dy < chN*^ /TogN llu II Bw,ll. 

y 2,n“ 


Proof: Let p be any real number, 2 < p < +“, and p be such that 1/p + 

* 2 

1/p = 1. Since u belongs to H (Q) , the trace 3u/3n belongs to 

1 / 2 

H ' (y); hence Lemma 3.5 gives 


B3u/3nB A < c/pBuB 
0,p,Y- 


2,fl 


Next , let us compute 


/ (3u/3n)(0,y)(w N - w h >(0,y)dy £ II 3 u/ 3 n II Q p ^Bw h - w N 


0,P ,Y 
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c/plul _|W “ W I ^ 

2, ft 0,p ,Y 


But, due to the definition (2.10) of , wj, is equal to -t^Wfi on y, so 

that, by (3.2) 


/ (3u/3n)(0,y)(w - w )(0,y)dy _< c/p hlul _lw I * . 

Y 2,fi l,p ,y 


Applying Corollary 3.4 gives 


/ (3u/3n)(0,y)(w - w )(0,y)dy_< c /p hN 

Y 


3 /2-1/p 


Hul _I w II , 

2,0 A/ ,T 


< c /p hN 1/2 N 1/p lul _lw . 

2,0 N 1/Z ’ 7 


Choosing p - logN, we obtain the desired result. 

Remark 3. 7 . Of course, the estimate (3.7) is not what we would wish, since 
the convergence is obtained only if the discretization parameters are linked 
by the following condition: 


(3.8) 


limthN 1/2 /T3iin = 0, 


(in fact, in (3.7) and in this condition, h can be replaced by h which is 
the greatest of the lengths of the edges of triangles K in T ^ contained 
in Y) . 


Remark 3.8: The estimate (3.7) is independent of k; indeed, we do not know 

how to improve it for large values of k. 
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3.2 The case of integral matching condition 

This case turns out to be simpler. We denote by the Hilbert pro- 

jection operator from \r (y ) onto x^. We have for any v in Hq(y) 


|V ‘ 1I h Vl 0,Y — ' V ' V'o.Y’ 


so that, for any v in H*(y) P i H*(Y) , 1 < t < k+1, 


(3.9) 


,V ' V'0,Yi Ch " V 't,Y- 


By interpolation, this inequality also holds for any v in 

Hq(y), 1/2 < l < 1. Finally, recalling that the interpolation space with 
index 1/2 between H*(y) and L 2 (y) is H^ 2 (y) (see [LM, Chap. 1, Th. 

1 /O 

11.7]) and denoting by H . N ^ the norm of Hqq (y), we also obtain 

1/2 1/2 ’ Y 
for any v in Hqq (y) 


(3.10) 


II v - ir.vll- < ch^ 2 Hvll * 

h °* T - 1/2 ,Y 


Now, we prove the following. 


Proposition 3,9 : For any function u in such that the func- 
tion u“ belongs to where l is a real number, 2 < l < k + 5/2, 

the following estimate holds for any w^ = (w h ,w N ) in 


(3.11) 


/ (3u/9n)(0,y)(w - w,)(0,y)dy < ch* 1 Hu H II w* I 

y 
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<H> i 

Proof: Let w c = (w,,w_ T ) be any element in vi; due to the definition 

o n N 6 

(2.11) of V^, coincides with tt^Wj^ on Y* We com P ute 

/ (3u/3n)(0,y)(w N - w h )(0,y)dy = / ( 3u/3n) (0,y)(w N - T h w N ) (0,y)dy 
Y Y 

= / [ ( 3 u/ 3 n) - ir h (3u/3n)(0,y)](w N - ir h w N )(0,y)dy, 

Y 

so that 

/ (3u/3n) (0 ,y) (w - w )(0,y)dy < ll(3u/3n) - ir (3u/3n)H «w - ir.w » . 

^ N n — h U,Y N h N 0,Y 

If i - 1/2 is not an integer, we note that 8u/3n belongs to 

Hq (Y)nH £-3 / 2 (Y) and that, since Wjj vanishes on 3ft + \y, w^j^. belongs 

to Hqq 2 (y). Applying (3.9) or (3.10) to bound the first term and (3.10) to 
bound the second one , we obtain 

/ (3u/3n)(0,y)(w N - w h )(0,y)dy _< ch £ 1 Idu/dnl^^ ^Iw^ * , 

Y 1/2 ,Y 

(with H3u/3nll . replaced by B3u/3nH * in the case £ = 2), 

' ,Y 1/2 ,Y 

and the result follows. The case where £ - 1/2 is an integer follows by 
an interpolation argument. 

Remark 3.10; Clearly, the estimate (3.11) is much better than (3.6), since it 
is independent of N. In fact, the term / ( 3u/3n) (0,y) (w N -w^)(0,y)dy 
goes to 0 whenever the discretization parameter h decreases to 0. 
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4. ANALYSIS OF THE APPROXIMATION ERROR 

We begin by recalling some properties of the approximation by finite 
element functions and by polynomials in the two-dimensional case. 

First, since for each K in the set is P^CK )-unisolvent , 

there exists an interpolation operator from {v e C°(ft ); v = 0 on 

3ft \y) into such that, for any function v continuous on ft and 

vanishing on 3ft - \Y, T^v is the only element of which satisfies 

(4.1) Vac S h n(ft"UY), (T h v)(a) = v(a). 

Q — 

Moreover, if the function v belongs to H (ft ) for a real number 
£, 2 _< £ £ k + 1, the following interpolation error holds [C, Thm. 3.1.5] 

(4.2) II v — T.vll _ ch^ m llvU _, m = 0 or 1. 

m,ft I, ft 

Next, we state the following result which can be derived in the same way 
as in [M2, Thm. 3.2]. 

Lemma 4.1 ; Let p be a real number 1 such that p - 1/2 is not an 
integer. There exists a projection operator n p from the space 
(v e H p (ft + ); v = 0 on 3ft + \y) onto {v N e Q^ft*" ); v N = 0 on 3ft + \y) 

such that, if a function v vanishing on 3ft + \Y belongs to H°(ft + ) for 

a real number o p, the following error estimate holds 

nv - n p vii . < cN y-a nvn , , o < v < p. 

N y,ft + “ a,ft + “ “ 


(4.3) 
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Now, we are going to approximate a function u of Hq(£ 2) by a func- 
* 

tion Vg , with in V^. In fact, we shall set 


(4.4) 


v 6 = (T h U " + V n N U+)> 


where will be chosen in so that satisfies the matching 

condition. Of course, the choice of depends on this condition. 


4 . 1 The case of ponctual matching condition . 

We immediately prove the following result. 

Proposition 4.2 : For any function u in Hg(£2) such that the pair u* 

belongs to )xH a (& + ) , where 1 and a are real numbers, 

2 £ l <_ k + 1 and a > 2, there exists a pair in VP such that 

(4.5) Hu* - v*H _< c{h* -1 llu"ll _ + (h 0-1 + N 1-a ) Hu + H }. 

£,a“ cr, a 

Proof: For any function z defined on n + , let us denote by *z the 

function defined on £2 by 


V- (x,y) e £2 , z(x,y) = z(-x,y). 


Next , we take 
(4.4) with q^ 


2 _< p <. inf{k + I , o } , p - 1/2 i IN, and choose 

«_ —p 4 1 1 

equal to “T^(u "" n U Clearly, we have 


* , 
as in 

at any point 


a of 
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(T^u ) (a) + q h (a) - u (a) - u + (a) + (l^u + )(a) = (II^u + )(a), 


so that belongs to V^. Moreover, we write 


III - Vg II < llu - T^u II _+ llq^ll + llu + - n^u + ll + 

1 , Q 1 , ft 1 , ft 


< nu - T h u n _ + ii ( id - T h )(u + - njju + )n _ + 2«u + - i^u + n . 

i,n i,ft” i,£2 


Finally, applying (4.2) and Lemma 4.1, we deduce 

Hu* - v*ll c{h £- 1 llu - n + h P_1 llu + - n^u + H + N 1 -Cr ilu + « } 

i,sf p,£2 _ a,fl 

c{h^ X |lu II _ + h * 5 *N P a nu + ll + + N* a nu + ll + }. 

£,fi a,0 a,n 

Applying the convexity inequality a3 a p /p + 3^/q (a 0,3 _> 0,1/p + 1/q = 
1), gives the proposition. 

4.2. The case of Integral matching condition . 

As far as the approximation error is concerned, this case is less 

simple. First, we introduce the following lifting operator from x^, 

into *h : for any v^ in xy,, Ryyvy, is equal to v^ on y and vanishes 

in any point a of S^\y. This amounts to stating that 


% ' Ja£ 5 t v h (a>,, a|Y C V Vh ' 

n h 


(4.6) 
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2 

We recall that ir. stands for the Hilbert projection operator from L (y) 

n 

onto x^* Our purpose is to prove a stability result for the operator 

We use now a technique due to M. Crouzeix and V. Thomee [CT]. We first 

introduce the subset 5? of the points of £, which are a vertex of a 

n n 

triangle K in and the subspace of 

(4.7) = {v h £ Va e 5°, v h (a) = 0}; 

1 2 

we denote by x^ the orthogonal subspace to x^ in x^ with respect to 

2 . 2 
the scalar product of L (y ) (when k is equal to 1, the space x^ is 

simply {0}, and the space x^ coincides with x h ) . Let tt*, i = 1 or 

2 i 

2, be the orthogonal projection from L (y) onto x^; we have of course 

(4.8) ^h “ ^h + "h* 

2 

Hence, we are going to prove a stability result successively for R^tt^ 

and Vh- 

a 

First, we recall some classical notation. Let K be the "reference" 

triangle with vertices (1,0), (0,1), and (0,0). For any triangle K in T^> 

A 

there exists an affine mapping F^ which maps K onto K; let B^ be the 

Jacobian matrix of F K , and h k be the diameter of K. Moreover, if the 

triangle K has an edge E contained in "y, we assume that F^ maps the 

A 

edge E with vertices (0,0) and (1,0) onto E, and we denote by Bg the 

A 

Jacobian matrix of the restriction of F R to E. For any function q on 

A 

K, we set: q = qoF R . 
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We denote by £ h the set of triangles K in T h which have an edge 

contained in . For any K in we set 


(4.9) 


Y r = Span{q a ,ae? h Pi 3K}. 


Lemma 4.3 ; Let K be a triangle of and let E be the edge of K 

contained in Y. Any polynomial q in Yj^ satisfies 


(4.10) 


'< l l i K , < , 0.1,E /l <E < C * 


Proof: Due to assumption (2.2), it is well-known (see [C, Thm. 3.1.2 and 

3.1.3] or [B, Lemma 2.3]) that, for any function q in H*(K), 

lq « , „ < cM B” 1 II | (det BL)| 1/2 «qll A < c'BqB 
1,K “ K K 1,K 1 , K 


and 


whence 


Hq II 


H q II 


0.1.E ■ l (det , ;< c 'V<’ n , 


* y 

E 


= | (det B )| 1/2 »qS A > c'h^BqB 
°’ E E 0,1,E“ * 


0,E 


Iql Iql /Bqll 2 < clql A lql Ji q II 2 
1,R- 0,1, E U,t ” 1,K 0,1, E 0,1 


Finally, if q belongs to Y R , q belongs to the space Y spanned by the 

Lagrange polynomials of the points (j/k, 0), 0 _< j _< k. Clearly, on this 

A A A A 

last finite-dimensional space, the semi-norm: q Bqll * and q Bq fl 

0,1, E 0, E 

are norms and the three norms B.B *,B.B * and n . || * are equiva- 

1,K 0, E 0,1, E 

lent, completing the proof of the Lemma. 
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We need the following Gagliardo-Nirenberg inequality. 


Lemma 4.4: Let s be a real number, 0 s < 1/2. Let E be the edge of a 

triangle of T^, contained in y. Any function w in H^(E) satisfies 


(4.11) 


nwii 


o,»,e 


< c llw H 


1/2 

s,E 


llw II 


1/2 

1-s ,E* 


Proof: Using the previous notation, we set w = woF K , and applying the 

classical Gagliardo-Nirenberg inequality, we obtain 


B wU 

0,",E 


A 1/2 A 1/2 
< ciwi 'tiwi ' A . 

s,E l-s,E 


Next, applying [C, Thm. 3.1.2 and 3.1.3] or [B, Lemma 2.3] gives 


nwn a < c l (det B E>r 1 / 2 |,wll o,E < c ' h K 1 / 2 |,wll o,E’ 

0^ E 

* —1 /2 1 /2 

nwn A £ ciBgilCdet Bg)| nwn 1 E < c'l ^ 7 iwij E 


Interpolating between these two inequalities, we have 


Iwl A < chj 1/2 flwll 
s,E 


s,E 


and 


II wl A 

1-s ,E 


c'h^ /2 S HwB 


1-s ,E* 


so that 


. „ , A . . . s/2-1 /4.1/4-S/2. 1/2 A/2 

llw II n m _ - Iwl * < civ tV Iwl _ ^llwll, 

°’ ’ E 0,-,E ~ 


1 II W II , _ « 

s,E 1-s ,E 


We are now in a position to prove the following lemma. 
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Lemma 4,5: Let s be a real number, 0 < s < 1/2* The operator 

satisfies the following stability property: for any function w in 


H 


l-s 


(Y), 


Vh 


(4.12) 



< c II w II ^ II w II 

- s,Y 


1/2 

l-s,Y* 


1 _g 

Proof : Let w be any function in H (y)* If ECy Is the edge of a 

2 

triangle K of there exists a polynomial in which is equal to ir^w 

on E and vanishes on TYE; using this polynomial and the definition of 
2 

we have 


llir h Wll 0,E = l w(y)<V*)(y)dy < llwll 0 ,» E ,ir h W “o, 1,E' 


Now, applying Lemma 4.3, we obtain 


2 2 2 2 

IIR.fr.wll , < II tt. wll _ _/ Hit, w» . < c' HwB. 

n h 1,K — h 0,E h 0,1, E — 0 ,»,y 


Using Lemma 4.4, we derive 


2 . . „ . 1/2 „ . 1/2 
IIR.ir.wl, . < clwl _ MwH , 
n h l,k — s,E l-s,E 


Adding up this inequality on all triangles K of and using the Cauchy- 

Schwarz inequality, we obtain (4.12). 

Next, we consider the operator ir^. Following [CT] and using the 
theory of orthogonal polynomials, we can easily find a polynomial ij) in 


P k ([0,l]) such that 
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(4.13) 


1 

V<|> e P ( [0, 1 ] )/4>(0) = «^(1 ) = 0, / <Ky)«|>(y)dy = 0, 

K 0 


i|>(0) = 0 and <J>(1) = 1 • 


We denote by (0,y^), 1 _< i _< M, the points of with 

0 < y^ < ... < y m < 1. Setting y 0 = 0 and y M+ i = 1, we denote by E i+1 / 2 > 

0 <_ i <_M, the edge with vertices (0,y i ) and (0,y i+1 ). Finally, we define 
the function , 1 _< i £ M, by 

4» ± (y) = <K(y - - y^)) if y e 

(4.14) ♦( (y - y 1 + i)/(yi - y 1+1 )) if y e ^ i+l/2 > 

0 elsewhere. 

It is an easy matter to see that the set (l]^, 1 i M) is a basis of xj. 
We next establish the following lemma. 


Lemma 4.6 : The following stability property holds for any i, 1 _< i < M, 

(4.15) IR.*.I < c. 

i,n~ 

Proof: As a result of (2.2), the support of is made of a finite number 

of triangles K. of T^, bounded independently of h (see [B, Rem. 3.2]). 
Let K be such a triangle, and assume that maps the point (0,1) onto 

(0,yi>« 
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1) If K meets Y only in (0,y^), ^h^i|k coincides with the Lagrange 
polynomial q associated with this point, hence 

» R h^i«i K < cl^ 1 1 1 (det B K ) | 1 /2 II q II „ < c'. 

1,K 


2) 'If K has an edge E contained in Y, we introduce the polynomial 

A 

of P^OO which coincides with on [0,l]x{0} and vanishes at any 

A 

point of S* which does not belong to E, i.e., at any point with coordi- 
K 

nates (j/k, i/k) , 0 j k, 1 k - j; then, coincides with 

* — 1 

i oF^ and 


HR. 




h^ii K < cl^lKdet B K )| 1/2 lhMI . < c'. 
* 1,K 


M 


completing the proof. 

2 1 11 

Let w be any function in L^(y). Setting ir,w = £ X . ij). , we note 

i=l 1 1 

that the vector of coefficients X = is given by 


(4.16) 


GX - y, 


where G = ^®ijh<i j<M i® t ' ie sc I uare matrix of order M defined by 


1 

(4.17) V(i,j), 1 < i,j < M, g y = / * 1 (y)* j (y)dy, 


and where y = is t * le vector ^ given by 


1 

Vi, 1 <_ i 1 M, y = / w(y)i|> (y)dy. 

0 


(4.18) 
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It is easy to see that the matrix G is symmetric and tridiagonal. Its 
coefficients have been computed in [CT, Lemma 2]: 


g ii = (y 1+1 “ + 2) and 


(4.19) 


g i>1+1 ■ (~l) k 1 (y 1+1 " y jL )/k(k + 1 )(k + 2). 


Using the arguments of [CT], we prove the following lemma. 


Lemma 4.7: With the notation (4.16) to (4.18), the following estimate holds 


(4.20) 




Proof: Let D be the diagonal matrix with coefficients d^ = g££, 

1 <_ i H. We have G = D(Id + K) , where K is a tridiagonal matrix with a 

zero diagonal; moreover, all the coefficients of K are bounded by l/(k + 
1). If II. II j denotes the matrix norm associated with the euclidean norm 

on if*, it has been computed in [CT] that, for any integer l > 1 , 


BK A II 2 < (21 + l) 1/2 (k + l) - *. 

Hence, we deduce 


« (Id + k) _1 ii 2 < 1 + iik a h 2 < 1 + (21 + l) 1/2 (k + i)' £ , 


so that 


II (Id + k) -1 ii 2 


is bounded by a constant depending only on k. This 


proves the lemma, since 
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X = (Id + K) _1 D -1 y. 


We are finally in a position to prove the following lemma. 


Lemma 4.8: Let s be a real number, 0 _< s < 1/2. The operator R^n^ 

satisfies the following stability property: for any function w in 

h 1 -s (y), 


(4.21) 


hrVwh < cHwi^iwhJ / 2 . 
1,0“ “ S>Y 1_S>Y 


Proof : 
have 


TT h W 


M 

I Xj*!. 


i=l 


Let w be any function in s (y). Setting 


1 M 

Vh" ' ,1, Wi 1 

1=1 


we note that each <|>^, 1 _< i M, is orthogonal in H^(fi~) to any 

'('j* 1 < .1 <M, ± but at most two, whence 


M 

IIR.ikwII 2 < 3 Y X 2 |R.f,| 2 
h h 1,0“ “ i=l 1 ^ 1 i,0" 


Lemma 4.6 then yields 


12 M 2 
IIR.tt^wII < c l X, . 

hh I,!)’" i-1 1 


we 


Using Lemma 4.7 with the notation (4.18) we obtain 


"Vh "" 2 - < c l wi/<y 1+ i _ y i> 2 - 

1,0 i=l 1 
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Noting that the support of ^ is E i-i /2 ^ E i+1 /2* we com P ute 


y. < ilwH_ , \ l|l Mn 1 (y,.i - M t \ B B o 1 rn 

i - 1 °» 1 » Y i+1 i_1 0 *“»(yi_i»y 1 + i> o,i,(o,d 


so that 


hr.iJwi 2 < c 

hh i,a"“ 


M 


M 


I » wl 'o co (v v 3 < c ' l " wll 0 - E 
i=l °» ’ Cy i-I’ y i+1 ; i=l °» ’ 1-1 


/2 


Thanks to Lemma 4*4, we derive 


1 2 M 

HR, TT.WlI < c y Hwll „ llwll , 

hh , A s >*i-i/2 s ’ E i— 1 /2 


1,0 


i=l 


and using the Cauchy -Schwarz inequality gives the result 


1 * 
Proposition 4.9 : For any function u in Hq( 0) such that the pair u 

belongs to * H CT (Q + ), where l and a are real numbers, 

* 1 

2 £. * <_ k + 1 and a 2, there exists a pair in such that 


(4.22) 


Hu* - vtl < c{h* -1 Hu~H + 

N 1 a Hu + H + }. 



£,0 

a,n 

Proof : 

Let 

s be a real number, 0 < s < 

1/2; we take p = 3/2-s. We 

choose 

* 

v 6 

as in (4*4) with q^ equal to 

Vh^ u+ " ^h u ~^* 711611 » we 

have 





W + q h ) • Vn u+ = °> 


so that belongs to . Moreover, we compute 
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lu* - v*n < lu - T.u i + ir. ir. (n£u + - T.u")n + iiu + - jl?u + h 

6 - h i,o" N h i,o“ * l ,o + 

< Hu - Tu II + HR, it, (u - T.u ) H + HR, tt (u + - IL?u + ) II 

h 1,0" hh h i,o" * 1,0 


+ p + 

+ lu - n„u h 

N 1,0 + 


Lemmas 4.5 and 4.8, together with (4.8), imply 


* * — — — — 1 /2 — —I / 2 

lu - v~l < lu - T.u II + Hu - T.u H llu - Tu H . 

6— h i ^ hs,Y h l-s,Y 


+ lu + - <uV^nu + - n N U+ Hl-s,Y + “ U+ " K U+ \ Q +* 


Using the trace theorem, we obtain 


* * — — — I /2 — # — 1/2 

lu — v.l < lu - T u H _ + llu - T,u H Hu — -Lu II 

0 “ n 1,0 s+1/2,0 3 /2-s , 0~ 


+ iu + - nJJu + H 1/2 lu + - nJ|u + H 1/2 + Hu + - II^u + H 

N s+l/2,fi + * 3/2-s,a + * 1,G + 


Finally, estimates (3.2) and (4.2), together with an interpolation argument, 
the trace theorem and Lemma 4.1 prove the result. 


Remark 4.10; Here also, the error is better for the integral matching condi- 
tion than for the ponctual one. Indeed, in (4.22), the two discretization 
parameters are involved in a completely independent way. 
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5. FINAL ESTIMATES AND CONCLUSION 

First, we recall an estimate which follows at once from a standard result 
in spectral methods [CQ2, Lemma 3.2] [MQ, Formula (3.22)]. 

2 + 

Lemma 5.1 : For any function f in L (ft) such that the function f 

belongs to H p (ft + ), where p is a real number > 1, the following estimate 

holds for any w = (w^, Wjg) in Vg 

(5.1) (f*,w*) - (f*,wj) 6 < cN 1 “ P Bf + l + #w N H + . 

c p,ft 0,ft 

Our main results are stated in the two following theorems. 

Theorem 5.2 ; Assume that the solution u of problem (1.1) is such that the 

pair u* belongs to H^(ft )xH°(fl + ), where o is a real number _> 2. 

Assume moreover that the function f of L (ft) is such that the function 
f"*" belongs to H p (ft + ), where p is a real number > 1. Then, in the case 
of the ponctual matching condition, the solutions u and u^ of problems 

(1.1) and (2.15) satisfy 

(5.2) 

Hu* - u*» < c{hN 1/2 /IogN Hu~« _ + (h 0-1 + N 1_CT ) lu + 1 + N 1_P lf + l }. 

2, ft a, ft p,ft 

Theorem 5.3 : Assume that the solution u of problem (1.1) is such that the 

pair u* belongs to H^(ft”)xH°(ft + ) , where f. and o are real numbers, 

2 < k + 1 and a > 2. Assume moreover that the function f of L (ft) 
is such that the function f* belongs to H p (ft + ), where p is a real 



I 
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number > 1. Then, in the case of the integral matching condition, the solu- 
tions u and Ug of problems (1.1) and (2.15) satisfy 

(5.3) Hu* - u^l _< c{h £_1 8u~« _ + N 1 II u + II + + N 1-p «f + ll + }. 

a, ft p,fi 

Proof: We set $ = (h,(N - l) -1 ). Of course, we apply Proposition 2.5 and, 

in (2.20), we choose v g = (\ ,v N ) equal to the pair defined in Proposi- 
tions 4.2 and 4.9 respectively, but with 6 replaced by 6. Since v N 

belongs to Q N _^(Ti + ) and the quadrature formula (2.12) is exact on all 

* 

polynomials of degree £ 2N-1, this implies that, for any w^ in V^, 

a (v *,w*) = a 6 (v*,w*). 


Then the estimates (5.2) and (5.3) follow from (2.20), Propositions 3.6 and 
3.9 respectively, Propositions 4.2 and 4.9 respectively, and Lemma 5.1. 

By a classical duality method, it is possible to derive an improved esti- 
mate for 8u - u-J rt . in the case of the integral matching condition, 
o 0 ,ft 

Proposition 5.4: Under the assumptions of Theorem 5.3, in the case of the 

integral matching condition, the solutions u and u^ of problems (1.1) 
and (2.15) satisfy 


(5.4) 


»u - u.l n < c{h* _1 (h + N *)lu « + N 1 °(h + N - 1 )lu + l , 

60 ’ fi “ A, ft" a,ft + 


+ N 1-p 8f + » + }. 

p , Q 




l 
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Proof : We have 


" u " Vo.o 


sup 2 / (u - u g )(x)g(x)dx/Hgll 0 

geL (ft) ft 


9 1 

Let g be any function in L (ft) . The unique solution w in H Q (ft) of 
the problem 


-Aw = g in ft. 


(5.5) 


w = 0 on 3ft, 


satisfies 


(5.6) Bwll 2,ft — C,gll 0,ft’ 

«n> 

Setting u g = (u h ,u N ) , we compute 


/ (u - u.)(x)g(x)dx = a(u - Ug,w ) + / (3w/3n) (0 ,y) (u^ - u^)(0,y)dy. 


ft 


Hence, for any w^ in V^, with $ = (h,(N - 1) *), we have 

6 


^ ^ ^ ^ ^ ^ ^ 

/ (u - u.)(x)g(x)dx = a(u - u., w - w g ) + (f ,w fi ) - (f ,w g ) 6 

ft 


+ / (3w/3n)(0,y)(u N - u h )(0,y)dy. 
Y 


£ 

Choosing w^ 


as defined in Proposition 4.9 and using Lemma 5.1, we obtain 
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/ (u - u^)(x)g(x)dx < c{(h + N *)llu - u.ll + N* P llf + H + }Hwll- _ 

a p,a ’ 

(5.7) 

+ / (3w/3n)(0,y)(u N - u h >(0,y)dy. 


It remains to estimate this last term. But we note that u h|y * s et ^ ua ^ 
to x.u„, so that 

n N 


/ (3w/3n)(0,y)(u N - u h >(0,y)dy 
Y 

- / [ (3w/3n) - x h (3w/3n)](0,y)(u N - ir h u N )(0,y)dy 

Y 

■ / [(3w/3n) - x. (3w/3n)] (0,y) [(u - x^u) - (id - x h )(u - u N )](0,y)dy 

Y 

< H(3w/3n) - x h (3w/3n)ll Q ^(llu - x h u» 0 ^ + II (id - x h )(u - u^Hq ^) . 

Using (3.9) and (3.10) yields 

/ (3w/3n)(0,y)(u N - u h )(0,y)dy £ c h 1//2 »3w/3nll A 
Y 1 /2 ,Y 

f . 1 - 1/2 - . , 1 / 2 , + „ , 

(h II u H + h lu-Uj.il ), 

z,n~ N 1,0 

which together with (5.6) and (5.7), gives (5.4). 

The detailed analysis we have performed allows us to compare the two 
algorithms, corresponding to different matching conditions. Indeed, whatever 
the regularity of the exact solution is, we obtain better convergence results 
in the case of the integral matching condition. Since we have already noted 
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that the computational cost of the two methods is of the same order, we think 
that this last algorithm has to be preferred. Numerical tests^ [KP] which are 
currently being implemented should confirm the theoretical results. 

As already stated in this paper, we are only concerned with a model prob- 
lem on a model domain. However, in this very simple example, it turns out 

that the order of accuracy in the finite element domain is simultaneously re- 
stricted by the degree of polynomials and by the regularity of the solution, 
while in the spectral domain it is only limited by the regularity of the solu- 
tion. That is why we believe that, in more general problems, the finite ele- 
ment domain must be chosen in such a way that it contains a neighborhood of 
both the singularities of the solution (in the case of hyperbolic equations 
with shock waves for instance) and the singularities of the boundary of the 
domain (for instance, corners of polygons which induce singularities of the 
solution even if the righthand member is very smooth). Then, local refine- 
ments of the mesh can be applied to improve the convergence, in a much simpler 

way than for the p-version of finite elements. These techniques are presently 

3 

being developed by the second author. 
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